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ABSTRACT 

We have measured the intrinsic disk cllipticity for 7 nearby, nearly face-on spiral galaxies by combining 
Densepak integral-field spectroscopy with /-band imaging from the WIYN telescope. Initially assuming 
an axisymmetric model, we determine kinematic inclinations and position angles from Ha velocity fields, 
and photometric axis ratios and position angles from imaging data. We interpret the observed disparities 
between kinematic and photometric disk parameters in terms of an intrinsic non-zero ellipticity, e. The 
mean ellipticity of our sample is 0.05. If the majority of disk galaxies have such intrinsic axis ratios, this 
would account for roughly 50% of the scatter in the Tully-Fisher relation. This result, in turn, places 
tighter constraints on other sources of scatter in this relation, the most astrophysically compelling of 
which is galaxy mass-to-light ratios. 

Subject headings: galaxies: structure - galaxies: kinematics 



1. INTRODUCTION 

A large fraction of galaxy disks appears to be non- 
axisymmetric. Photometric studies of face-on galaxies 
find the light distri bution in 30-50% of spirals is lopsided 
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Zaritoky fe Rix, 1007 j- Rudnick fe R i* 



1998 )-Q¥- has higher order distortions (m > 2; | Rix fc Zarit 



sky, 1995| ; |Kornreich et al., 1998| ; |Conselice et al 2()()0| ) 
A similar fraction of galaxie s haslopsided HI distributions 
or kinematic asymmetries ( Baldwin et al., 1980 [_ Bosma 



1981} iRichter fe Sancisi, 1994| ; |Hayncs et al., 1998| ; ^wa- 
ters ct al, 1999| ). 

Disk cllipticity (to = 2) may also be common. Few 
disk galaxies appear round on the sky; the distribution 
of apparent axis ratios is best modeled by ran domly ori- 
ented disks with intr i nsic ellipticities e ~ 0.1 (|Binney fe 



de Vaucouleurs, 1981 ; Gr0sbol, 1985 ; Huizinga fe van Al-_ 
bada, 1992j ; jLambas et al., 1992| ). These statistical studies 



constrain mean ellipticity, but do not probe the intrinsic 
ellipticity of individual galaxies. 

If galaxy disks are elongated, errors in photometric es- 
timates of inclinations and position angles (PAs) will be 
introduced, thereby contributing to the scatter in the re- 
l ationship between ga la xy luminosity and rotatio n speed 
dTully fe Fisher, 1977[ ). [Franx fe de Zeeuw (199~2] ) used a 
specific, non-rotating velocity field model to show that if 
e = 0.1 on average, systematic errors could account for all 
the observed scatter in the Tully-Fisher (hereafter, TF) 
relation. However, the true contribution to this scatter 
depends on the distribution of intrinsic axis ratios. Iden- 
tifying the individual ellipticities of galaxies would limit 
other sources of astrophysical scatter in the TF relation 
and reduce scatter. 

Kinematic maps can be used to estimate ellipticity for 
individual galaxy disks. One approach involves fitting such 
maps with elliptic orbits, from which a quantity related to 
the ellipticity of the potential can be estimated: e po t sin 20. 
This quantity varies from 0.001 - 0.07 for a sample of 9 



galaxies (Schoenmakers et al., 1997; Schoenmakers, 1999). 
Unfortunately, the phase angle <fi is not directly measur- 
able. Such studies have been rare because they require 
observationally expensive, hig h resolution, high signal-to 
noise velocity fiel d maps (e.g. Teuben et al., 1986; Korn 



reich et al., 2000) 



Here we develop a new, efficient method of estimating 
intrinsic galaxy disk ellipticities using optical kinematic 
maps and photometric indices. Determining unique solu- 
tions for e and <p requires measurements of both kinematic 
and photometric inclination and PAs for nearly face-on 
spiral disks. When the disk major axis is rotated by 
a phase angle <p away from the line of nodes, the kine- 
matic and photometric axes will appear misaligned only 
if cosi < 1 — e. If e « 0.1, as suggested by earlier 
studies, substantial misalignment of kinematic and pho- 
tometric PAs (> 30°) will occur only for galaxies with 
i < 30°. To demonstrate our method we present the de- 
tailed analysis of WIYN integral field Ha velocity fields 
and / band images for UGC 4380, and measurements of e 
for an additional 6, nearly face-on galaxies based on the 
same method. 



2. SAMPLE AND OBSERVATIONS 

Our sample of seven galaxies is a subset of 69 selected 
from the Principal Galaxy Catalog (PGC; Paturcl ct al.J 
1997) in regions of low Galactic extinction with normal 
surface-brightnesses, axis ratios corresponding to < 30° 
inclination for a circular disk, diameters D25 between 0.75 
and 1.5 arcminutes (to facilitate spectros copic follow-up) 



and H ubble types between Sb and Sc (see Andersen et al 
2001). We excluded galaxies with obvious bars, clearly 



visible non-axisymmetric structure, and foreground stars 
within 2-4 disk scale lengths using the Second Palomar 
Observatory Sky Survey. 

We used the WIYN S2KB imager (a 2048 2 CCD with 
0.195 arcsec/pixel) to obtain deep J-band images of our 
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7 targets between May 10-14, 1999 with seeing of ~1.1 
arcseconds FWHM, or 10% of 1 scale length. We re- 
duced these images using standard procedures, yielding 
sufficient signal-to-noise to establish galaxy isophotes to 
~ 4.5 scale lengths. We used flux-calibrated, i?-band mag- 
nitudes obtained with the IGI/TK4 imager on the McDon- 
ald Observatory 107" telescope to check that these galaxies 
had close to "Freeman" disks (S-band central disk surface 
brightness of 21.7 mag arcsec -2 ). For example, assuming 
B-R = 1.15 for the Sc-type galaxy UGC 4380, we find 
jUo(-B) = 21.5 mag arcsec" 2 . 

Spectral images of our 7 targets were obtained using 
Densepak on the WIYN telescope from portions of runs on 
May 20-21, 1998, January 20-22, 1999, and March 27-28, 
1999. Densepak is a fiber-optic array used for integral field 
spectroscopy containing 86 active, 3-arcsecond fibers (sep- 
arated by 4 arcseconds) arranged in a 7 x 13 staggered grid 
subtending an area of 30 x 45 arcseconds ( Bardcn et al.J 
1998[j] Four additional "sky" fibers are spaced around this 
grid roughly one arcminute from the grid center. Densepak 
feeds the WIYN Bench Spectrograph, which we used with 
an echelle grating to cover 6600A < AA < 7000A with a 
dispersion of 0.19 A per pixel and a FWHM spectral reso- 
lution of 0.51 A. Basic spectral reductions were done using 
the NOAO IRAF package dohydra. Thorium- Argon lamps 
were used for wavelength calibration. 

Since galaxy rotation curves peak at roughly t wo pho 



tometric scale lengths (Courteau & Rlx, 1999; Willick 
19990 we used Densepak to map out to ~ 3 scale lengths 
per galaxy (typically ~ 30 arcsec). For UGC 4380, which 
has a /-band scale length of 9.7 arcseconds, two Densepak 
pointings were required to cover a square of roughly 45 x45 
square arcseconds. Two 30-minute exposures were made 
at each Densepak position to facilitate cosmic ray removal. 
The Ha line flux was typically measured at a signal to noise 
of ~ 10 at the edge of our observed field. 

3. MEASUREMENTS 

To measure photometric axis ratios and PAs we use 
the STSDAS ISOPHOTE package ellipse on images with 
stars masked. Figure 1 shows the radial dependence of 
PA, apparent ellipticity (1-b/a), and an azimuthal vari- 
ance statistic for UGC 4380. Inside 29 arcseconds, the 
variance statistic is greater than unity showing that spiral 
arm structure affects the ellipse fits. Hence we measure 
the photometric axis ratio and PA at radii greater than 
29 arcseconds, where the axis ratio and PA are constant, 
consistent with our estimate that spiral structure has a 
minimal impact. 

To determine kinematic inclination and PAs, we first 
produce Ha emission line maps from the 86 field-flattened 
and sky-subtracted spectra from each pointing. We mea- 
sure velocity centroids by fitting a Gaussian line profile to 
each Ha line, and assign a spatial position based on the 
fiber geometry of Densepak and the offsets used to make 
multiple pointings. Figure 2 shows a polynomial surface 
fit to these discrete velocity measurements for UGC 4380. 
UGC 4380 clearly is an inclined galaxy with a center co- 
incident with the photometric center. The formally well- 
determined PA and peak rotation velocity (Table 1) agree 
with visual inspection of the data (Figure 2). Single dish 
HI observa tions of UGC 4380 list a IU50 velocity widths of 
118 km/s (Haynes et al., 198§|). The W50 width measured 



from our Densepak spectra is also 118 km/s, implying we 
observed the peak of the rotation curve. 

To extract the PA, peak velocity and inclination, we fit 
the velocities expected at each fiber position using a sim- 
ple model consisting of concentric and coplanar circular 
orbits with V(r) = V c tanh(r//i), where the terminal cir- 
cular velocity V c and h are free parameters. The other free 
variables are inclination, PA, center, and central velocity. 
Our velocity fields exhibit only small deviations from this 
simple model, as illustrated for UGC 4380 in the middle 
panel of Figure 2. The best fitting model was determined 
from x 2 minimization (downhill simplex method) based 
on comparing the measured velocity centroids, fiber by 
fiber with the smooth model velocity field sampled with 
the Densepak footprint. The standard deviation in the 
fit residuals is 4.3 km/s. More elaborate radial velocity 
functions do not provide smaller residuals yet have more 
independent variables. Formal confidence limits (CL) were 
placed on these quantities by determining surfaces of con- 
stant x 2 ■ Table 1 contains our measurements of axis ra- 
tio, photometric PA and kinematic inclination and PA. In 
UGC 4380, the photometric and kinematic PA differ by 
9.5° ± 3.9° (68% CL). This disk is unlikely to be intrinsi- 
cally circular. 

4. ANALYSIS AND RESULTS 

We can estimate the ellipticity required to produce 
the above level of discrepancy by modeling the observed 
galaxy as a disk with intrinsic, photometric ellipticity e 
(right panel, Figure 2) but circular orbits (justified be- 
low) . The vector describing the ellipse in the galaxy plane 
is f ' = [—(1 — e) sin 9, cos 9] , where 9 = is the major 
axis. To project the ellipse onto the (observed) sky plane, 
it is inclined by an angle i about the line of nodes at an 
angle <p from the intrinsic major axis. In the plane of the 
sky, the ellipse describing the galaxy isophotes is given by 
a transformation matrix inv olving 9, i, e, and <j>. 

Franx & de Zeeuw (ryi?) found that for orbits in a flat 
but elliptic disk, fitting the velocity field with tilted circu- 
lar rings yields the correct disk orientation to first order 
in ellipticity. Accordingly, we take the kinematic inclina- 
tion and PAs to represent the true inclination and PA of 
the disk; e and <f> can then be determined given measure- 
ments of the apparent photometric axis ratio (b/a), the 
kinematic inclination (i), and the difference between pho- 
tometric and kinematic PAs (ip). 

When the true inclination i is closer to face-on than the 
photometric inclination, as is the case for UGC 4380, the 
galaxy must be at least as flattened as e > 1 — |cos« _1 . 
The intrinsic flattening must be even larger if the true 
major axis does not lie in the plane of the sky. This lower 
limit for UGC 4380, e > 0.06 ± 0.03, is inconsistent with 
a purely circular disk. By utilizing the 9.5° ± 3.9° mis- 
alignment of the photometric and kinematic PAs, a better 
estimate of the ellipticity can be obtained. 

In general, three equations relate the three observables 
b/a, i, and ip to the three unknowns e, <f>, and the angle 9 
at the apparent major axis. Using the Newton-Raphson 
method for nonlinear equations, we find e = 0.07t8;{j[} 
(99% CL) for UGC 4380, a solution inconsistent with a 
circular disk. Derived values for the seven galaxies, plot- 
ted in Figure 4 and listed in Table 1, range from e = 0.02 to 
0.20. Three galaxies (UGC 4380, NGC 2794, UGC 5274) 
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are inconsistent with having circular disks at the 99% CL; 
two galaxies (NGC 3890, NGC 5123) have ellipticities in- 
consistent with circular disks at the 95% CL. Only two 
galaxies (UGC 7208 and UGC 10436) are consistent with 
having circular disks within their 68% CL. The galaxies 
with the highest derived ellipticity e have near 90°; the 
line of nodes is almost perpendicular to the true major 
axis. This is consistent with our selection of round, ap- 
parently face-on systems. The two galaxies with e > 0.1, 
NGC 2794 and UGC 5274, both show evidence for faint, 
interacting companions; NGC 2794 has an AGN. 

Finally, we have checked that potential sources of sys- 
tematic errors - lopsidedness, spiral structure, and warps - 
do not affect our ellipticity measurements appreciably. Ta- 
ble 1 contains the (A\) amplitude of the m — 1 component 
of the Fourier expansi on of the light profile as defined by 
flZaritsky fc Rix, 1997j ). With the exception of NGC 5123 
these galaxies arc not significantly lopsided ((j4i) < 0.2), 
nor d oes (Aj ) correlate with our derived values for tjj or e. 



None of these galaxies show any sign of kinematic asym- 
metry. A 180°-rotational asymmetry measure, defined as 
^180 = I y,(Vobs + Vi8o)l/|2 Y. V bs\, akin to the ph otomet- 
ric asymmetry parameter of |Conselice et al. (2000D , yields 
close to null values for all. 

Strong spiral structure drives photometric PAs to 
change with radius R at a rate of (9PA/<91ogi? = cot6> p , 
where 9 p is the pitch angle of the arms. A warp in the disk 
also would be manifest as a twisting PA with radius. As 
discu ssed in §3, we make our photometric measurements 



between 3-4 scale- lengths where we find the photometric 
PA and axis ratio are constant, and our azimuthal vari- 
ance statistic (Figure 1) corroborates that spiral structure 
is no longer a dominant photometric effect. The velocity 
fields exhibit no residual structure correlating with radius 
or azimuth, despite the fact that these measurements are 
within 3 scale-lengths (in contrast to our photometric mea- 
surements). We find no evidence for twisting PAs in the 
velocity fields which is not surprising since even galaxies 
with strong outer warps usu ally have pla nar HI distribu- 
tions within 3 scale lengths ( Briggs, 1990Q . 



We have demonstrated that high-quality Ha velocity 
maps and I band images of nearly face-on galaxies can 
be used to exclude the hypothesis that galaxy disks are in- 
trinsically free ofm = 2 (elliptic) distortions. Our method 
for estimating the deviation from circularity suggests UGC 
4380 has an intrinsic ellipticity of ~7%. This is one of the 
first unambiguous detections of disk ellipticity for an in- 
dividual galaxy. If disks are intrinsically elliptic, the pho- 
tometric and kinematic axes will be misaligned in general, 
and the inclinations derived from isophote shapes will dif- 
fer from the kinematic inclinations. Deviations are partic- 
ularly large for face-on galaxies. WIYN/Densepak echclle 
spectroscopy and optical imaging are efficient means to 
estimate disk ellipticity for large samples of such systems. 

The intrinsic ellipticity for our 7 galaxies is e = 0.05 ± 
0.01 (error weighted mean and uncertainty), inconsistent 
with purely circular disks. Most of our targets are nor- 
mal, intermediate-type spirals, typic al of those selected 
for TF surveys. According to the Franx fc de Zccuw| 
(19921 ), e = 0.05 should produce -50% of the observed 
TF scatter in red and near-infrared bands. The cause of 
disk ellipticity is currently unclear, e.g., halo triaxiality 
or non-uniform matter accretion could be responsible. If 
larger samples show disk ellipticity at these levels, this 
implies variations in either disk mass fractions or mass-to- 
light ratio of today's spiral galaxies must contribute under 
0.1-0.2 mag dispersion in the luminosity of galaxies at a 
given rotati on speed. High resolution co s mological simu- 
lation s (e.g. Steinmctz & Navarro, 199G; van den BoschJ 



2000) indicate that variance in disk mass fraction should 



5. SUMMARY AND DISCUSSION 



only induce scatter along the TF relation, because halo 
contraction is greater for larger disk masses. Assuming 
the remaining scatter in the observed TF relation is dom- 
inated by differences in spiral galaxy mass-to-light ratios, 
this modest variance places strong constraints on the for- 
mation histories of such systems. 



We thank S. Barden at NO AO and D. Sawyer at WIYN 
for making Densepak happen. This research was sup- 
ported by the UW-Madison Graduate School, and NSF 
grants AST-9618849 and AST-9970780 (MAB). 
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Table 1 



ID 


b/a 


PAphot 


ikin 


PA kin 


4> 




e 






(dog) 


(dog) 


(dog) 


(dog) 


(dog) 




UGC 4380 


0.90±0.02 


43±4 


16.1±2.3 


33.5±1.8 


10±9 


15±15 


0.07 


NGC 2794 


0.85±0.04 


408±3 


20.5±1.4 


327.9±0.9 


80±3 


83±4 


0.20 


UGC 5274 


0.94±0.02 


78±7 


32.2±3.8 


9.5±1.4 


69±7 


84±5 


0.20 


NGC 3890 


0.96±0.02 


264±16 


23.5±1.8 


225.6±1.4 


38±16 


73±13 


0.08 


UGC 7208 


0.94±0.01 


325±9 


23.8±2.3 


330.4±1.2 


-5.4±9 


101±90 


0.03 


NGC 5123 


0.92±0.01 


163±4 


22.2±1.1 


153.0±0.5 


10±4 


44±27 


0.03 


UGC 10436 


0.87±0.03 


270±3 


30.6±1.4 


267.1±0.7 


3±3 


64 ±90 


0.02 



e confidence limits 



WT7- 



(Ai) 



+0.04 
-0.04 
+0.07 
-0.07 
+0.07 
-0.06 
+0.05 
-0.04 
+0.04 
-0.03 
+0.02 
-0.02 
+0.06 
-0.02 



+0.06 
-0.06 
+0.11 
-0.11 
+0.12 
-0.10 
+0.07 
-0.07 
+0.06 
-0.03 
+0.04 
-0.03 
+0.09 
-0.02 



+0.08 
-0.06 
+0.13 
-0.13 
+0.14 
-0.11 
+0.08 
-0.08 
+0.08 
-0.03 
+0.04 
-0.03 
+0.11 
-0.02 



0.06±0.02 
0.04±0.02 
0.07±0.02 
0.13±0.06 
0.12±0.05 
0.25±0.10 
0.09±0.05 
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Fig. 1. — The radial dependencies of photometric position angle (PA; top panel), ellipticity (1-b/a; middle panel) of UGC 4380's /-band 
isophotes. The bottom panel shows the azimuthal surface-brightness variance around each isophotal ellipse, normalized by the expected 
shot-noise (source noise, cr s , and sky plus read-noise, at,)- This normalized variance is large (> 1) where spiral structure contributes to the 
overall variance budget. We measure the PA and b/a between 29 and 39 arcseconds where both the PA and axis ratio remain constant, the 
azimuthal surface-brightness variance is consistent with shot noise, and the signal to noise per pixel is greater than one. 

Fig. 2. — The Dcnscpak Ha velocity field of UGC 4380 (left panel) and the residuals between the velocity field and a simple model (middle 
panel; see text). Both are smoothed and interpolated using a polynomial surface. Solid, heavy, and dashed lines are positive, zero, and 
negative velocities, respectively, relative to the model systemic velocity (left panel) or model velocity field (middle panel). The dash-dotted 
lines (left panel) represent the isophotal annulus determined from Figure 1 from which the photometric b/a and PA are derived (solid lines 
in annulus indicate photometric major and minor axes). The right, summary panel shows a schematic diagram of an isophote for UGC 4380 
assuming its disk has an intrinsic ellipticity of e = 0.07. The x-y plane has been rotated 33.5° from North to match the kinematic position 
angle of UGC 4380. The solid ellipse represents the true shape of the disk seen face-on; the major axis lies at an angle <j> = 15° from the line 
of nodes (the kinematic major, or y axis). The dashed ellipse represents the apparent shape after inclining the elliptic disk by i = 16.1°; it 
has an apparent axis ratio of b/a = 0.90. The kinematic and photometric PAs differ by if) = 9.5°. 

Fig. 3. — Solutions for ellipticity e and angle </> for seven sample galaxies. The 68% confidence limits (CL), shown as contours, are derived 
from the estimated measurement errors on b/a, tf) and i listed in Table 1. UGC 7208 (dark-dash) and UGC 10436 (light dash) are consistent 
with a circular disk within the 68% CL. 
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